Non-Magnetic Impurity induced in-gap bound states in two band s± superconductors 
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In this paper we study the effect of single non-magnetic impurity in two-band s-wave supercon- 
ductors where the two s-wave order parameters have relative phase S ^ between them. We show 
that in-gap bound states are always induced by non-magnetic impurities when 5 — tv (s±-wave 
superconductors). The bound state is a consequence of the topological nature of the corresponding 
Bogoliubov-de Gennes equation. 



With the discovery of the Iron-based (pnictides) su- 
perconductors, the superconductivity characterized by 
more than one order parameters, i.e. the multi-gap su- 
perconductors, becomes a hot topic. Band structure 
calculations indicate that the material has a quasi-two- 
dimensional electronic structure, with four bands cen- 
tered around the T- and M- points in the Brillouin zone 
contributing to the Fermi surface. It has been proposed 
that the superconducting order parameters in this multi- 
band materials have s-wave symmetry, but with opposite 
sign between bands centered at F- and M-pointsP,[3,[l]. 

The effect of impurities in this class of materials have 
been an issue of interests. NMr[J and lower critical field 
dataQ seems to indicate the existence of nodes in the 
superconducting order parameter while the observations 
in Angle-resolved photoelectron spectroscopy 0, 0] favor 
node-less gaps. One possible solution to this controversy 
is that the order-parameters are s-wave-like, but large 
number of in gap states are induced in the material be- 
cause of the frustrated sign of order-parameters between 
the bands. Indeed, such a scenario has received supports 
from self-consistent-Born type calculations where in gap 
states are found to appear easilvfsllgl.fioj. 

To have a more precise understanding of the physi- 
cal effect of non-magnetic impurities and the origin of 
in-gap states, we consider in this paper the effect of a 
single non-magnetic impurity on two-band s-wave super- 
conductors. The electrons in the two band are coupled 
through a Josephson effect which determines the rela- 
tive sign between the order parameters 14| . The effect 
of impurities is studied by analyzing the corresponding 
Bardeen-Schrieffer-Cooper (BCS) theory. In particular 
we want to understand (1) how the superconducting gap 
is affected by the impurity and (2) the conditions under 
which in-gapped states arc formed. We note that the 
problem of single impurity has also been tackled numeri- 
cally in particular tight-binding models used to simulate 
iron pnictides 11 , [isj . Our analytical work here is 



independent of the microscopic details of the models and 
provides results that are complementary to the above nu- 
merical works. 



formulation 

We start with the Path integral formulation of BCS 
theory. The system we consider is characterized by the 
BCS action, S = So + Sj, where 



i.k 
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A, 



Oi 
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1,2 is the band index and 



k — (k,iujn)- is the energy dispersion for elec- 

trons in band i and Ca , are spin-a electron annihi- 
lation(creation) operators. Sq is a sum of two bulk BCS 
mean-field actions describing two superconducting bands 
coupled only by Josephson interaction. Aq^ is the super- 
conducting gap when impurities are absent. The effect 
of a single non-magnetic impurity is represented by Si, 
where 



Si = 



A* — r/- 



(2) 



where il = volume of system and ^'i(i(jj„) = J^k^iik)- 
Si describes the effects of an impurity scattering po- 
tential t/(f) ~ S''{r)J2i^,J=l,2),aU^3cl^''^ir^c^'\r), where 
Uij 's are the scattering matrix element between bands i 
and J and A^ ^ = Sij A^ is the induced change in local su- 
perconducting gap as a result of the impurity scattering 
potential. We have approximated the induced change in 
gap to be of form Ai{r) ~ S'^{r)Ai here, consistent with 
our simplified form of impurity scattering potential. 

The superconducting order parameters are determined 
by the mean-field equation 



A,(f) = Ao. + A,(?^) 



where (cj"'' (r)cj"'^ (r)) is the pairing amplitude between 
electrons in j*'' band, Vu represents the pairing inter- 
action between electrons in band i and V12 = V21 is the 
Josepshon coupling between the pairing order parameters 
in the two bands. 



{c\'\r)cfif)), (3) 
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The fermion fields in S can be integrated out to obtain 
an effective action Sgff in terms of Uij and A^, we obtain 

Seff = lndct(il/o) + lndct(l + GoMi{U, A)), 

where Indet(Mo) is coming from the the mean-field BCS 
action in the absence of impurity and 



where 



goiiiujn) = 



-iuJn Aoi 



(4b) 

is the (on-site) Nambu matrix Green's function for band 
i electrons in the absence of the impurity, oju >> 
|Aoi I, I A02I is the cutoff energy for attractive interactions 
(~ Debye frequency for phonon superconductors) and 
Ni{0) ~ Ep^ is the band i density of states at the Fermi 
level. We have assumed Ep » ujd, ^oi{2), Uij, ujn, etc. 
to justify using a constant density of states in Eq. ([4]). 
We also have 



solved analytically in the limit y,|f/'|,|Ao| and co'jj = 
ttN{0)uju << 1, which is the case for weakly-coupled 
BCS superconductors. In this case it is straightforward 
to show that 

A^-Ao\U'\^ + 0{{U',V',uj'j,t), 

and the effect of impurity is to reduce the gap amplitude 
at the impurity site. Correspondingly a bound state is 
induced at the impurity site which is determined by the 
equation 1 + Go{u!)Mi = 0, or 



V|Ao|2-t^2(i + [;'|2 ^ |^/|2) ^ 2Re{A*A') = 0. 



We sec that a solution ui < |Ao| exists when A^A < 0. 
The solution has energy uj Ao(l - 2|A[||2|[/'|4) >> 
Ao- |A| in the limit \U'\,V' ,uj'jj « 1. The bound state 
solution is a direct consequence of local suppression of 
superconducting order-parameter by the impurity which 
creates a local "potential well" in the system. The bound 
state has energy uj > local gap magnitude = Aq — |A| and 
is therefore not a true "in-gap" state. 



two-band situation 



The free energy associated with the impurity is 

Fi = J2 lndet(l + GoA/i([/,A)) (6) 

+ Y,{A*.^ + A*){V-%{A„,+A,) 



and the mean-field equation for At can be obtained by 
minimizing the free energy with respect to A^. We obtain 



^ [(l + Go(ic^„)Afi)~'Go (iun)]- ■ = Ai + Aoi 



/9 



where (ji, J2) = {2j - 1, 2j) for j = 1, 2. 



(7) 



single-band case 

It is helpful to first consider the situation of single band 
superconductor. In this case the mean-field equation for 
A is 



TrN(0)(A' + - ^2 



A'o+A' 



V|AoP-(ic^„) = 



(8) 

where X' = TrN{0)X, where X = [/, y, A, Aq, wd. The 
BCS mean-field equation in the absence of impurity is 
recovered if we set U' — A' = 0. Eq. ^ can be 



Next we consider the two-band situation. To see the 
new physics associated with the appearance of multiple 
bands we first consider bound states assuming A^ = 0, 
i.e. there is no induced local change in the gap ampli- 
tudes. 

In this case we obtain after some algebra 



det(l + Go(w)Afi) = 1 + 



\Ui2\H\Ui2\'-^UiiUL) 
(l + |C/{iP)(l + |C/^2P) 



(9) 



2{lo^ -\Aoi\\Ao2\ cos S)\U{2 



(1 -f- \U[,m + |[/^2P)V(|AoiP-^^)(|Ao2P-^2) 



where [//^ = ir^jNi {0)Nj {0)Uij and S is the relative 
phase between the two order parameters Agi and A02. 
The solutions to the equation det(l -I- Go{uj)M2) ~ 
can be obtained easily since the equation is quadratic in 
a;2. We are interested at the bound state solution with 
u! < mm( I Aoi I, IA02I). Notice that the equation has no 
solution when U12 = 0, consistent with what we observe 
in the single-band case. 

Assuming that IA02I > |Aoi|, it is easy to see from 
Eq. that bound state solution with energy uj < |Aoi| 
exists only when |Aoi| > |Ao2|cos(S. In particular, no 
bound state solution exists when 6 = 0. However bound 
state exists when 6 = n, when the two superconducting 
order parameters are out of phase, even when there is no 
induced local changes in the order parameters. Solving 
the equation we obtain 

= 2(1 -4r^) + + 8r2|Aoi||Ao2| (10) 

-(|Aoi| + |Ao2|)yCA;;^ 



|Aoi|)2-M6r2|Aoi||Ao 
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where 



where 



(1 + |c/{iP)(i + m^P) - \ui,\H\u[,\^ - 2Ui,u^,) ■ 

We note that bound state sohitions always exist in the 
ergy 



hniit of small IC/i'tP (Born limit). The solution has en- 



-lA 



oil 



-2r- 



(|Aoi 



^02 I 



(|Ac 



lAoi 



■IIAoil +0(r4), 



in the limit |Ao2| - |Aoi| >> 4rA/|Aoi||Ao2| and be- 
comes uj — |Aoi| ^ — 2|Aoi|r in the opposite limit 
IA02I - |Aoi| « 4rv/|Aoi||Ao2|. More generally, it is 
straightforward to show that solutions with w > ex- 
ists when < 1/4 and no solution exists at > 1/4. 
Therefore there is an intermediate range of parameters 
U's where bound state solution does not exist. At around 
the resonance point = 1/4 — S we obtain 



UJ 



A6- 



AoiP + |Ao2|- 



Our result indicates that the existence of in-gap state 
and the corresponding bound state energy depends on 
the particular form of impurity scattering potential which 
determines the parameters Uii,U22 and Ui2- Assuming 
the f/'s are all proportional to each other we find that 
in the strong scattering (Unitary) limit iC^/jP — * 00, 
—^ and the bound state energy approaches zero 
asymptotically. This result is in qualitative agreement 
with tight-binding calculations 11 , [13, [31 where shal- 



low in-gap states are found to exist easily in two-band 
superconductors with order parameters of opposite sign. 
Our non-niodel-specific result suggests that uj bound 
states are in general allowed in two band superconductors 
with frustrated sign between order parameters. 

To examine whether the in-gap state is robust against 
changes in the superconducting order parameters we con- 
sider the case of symmetric bands with gQi{iuj) = (702 (*'^) 
and Aoi = Aoe**^*, i.e. the two bands differ only in the 
phase of the order parameters. To simplify the problem 
further we set J7ii = U22 = so that 



Ml 



Ae 



V 
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Ae'^'i 
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U \ 
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Ke-'^^ / 



(11) 



where A and 9i are to be solved self-consistently from the 
mean- field equation ([7]). The determinant det(l-l-GoMi) 
can still be computed analytically in this case. We obtain 
after lengthy algebra 

det(l -I- Go{iuj)Mi) = A{iuj) + S(iw)[cos(6ii - ipilU) 

+ C0S(6I2 - 02)] + C(jLj) [C0S(6I2 - 0l) + C0s(6'i - 02)] 

-h£i(za;) cos(0i - 02) + 2|C/{2p|AT cos(6li - 6I2) 
+E{iu}) cos(6'i - 0i) cos(6'2 - ^2), 



A{ioj) = {\ + \K'\f + \U'\^ + 2 



{iojf\U'\'' 



2AoAV|AoP-(zu;)2[i + |A'| 
(*^)2-|AoP 

2AoAV|AoP-(»c.)2|[/f 



(zu;)2-|Aop 



(*c.)2-|AoP 



where X' = ttN{Q)X as before. 

The mean-field equation can be solved in the weak im- 
purity scattering limit |C/'|, |Aq|, << 1. Keeping 
only terms to order |C/'P and A' in the mean- field equa- 
tion, we obtain 6^ ~ (f>i{i ~ 1,2) and 

A~ |C/f Aocos(0i-02). 

Notice that the superconducting order parameter is en- 
hanced by scattering between the two bands if 0i — 02 
but is suppressed by scattering if 0i — 02 = ±7r, suggest- 
ing that non-magnetic impurity induces a local ferromag- 
netic Josephson coupling between the superconducting 
order parameters which disfavors a s± state 0. 

We next examine the solution(s) to the equation 
det(l + Go(w)Mi) = with < |Aop. Defining 

= Aq—uj'^ we find that a zero at some < j/o < Aq im- 
plies an in-gap bound state. Solving the equation (with 

= 0i(* = 1; 2)) we obtain i/q = {—(5 -\- j)/a, where 

(3 = 2AoA'{l + |A'|2 + |C/'|2 cos(0i - 02), 

a = l+2(|C/f +|Af ) + [;7'4+|A'|4+2|C/'nAf cos(0i-02)], 
and 
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2 - 2|C/f A2[1 + |C/f + |A'|2]2[1 - cos(0i - 02) 



It is easy to see that 



yo 



2AoA' 



i + \u- 



l\2 



,/|2 



(13a) 



for 01 = 02 and a y > Q solution does not exist for small 
|?7{2p where AqA' > 0. The situation is very different 
for 01 — 02 = ±7r. In this case 



2/0 



2Ar 



|C/'|(|(1 + |A'|2 + |C/'|2)-A'(1 + |A'|2-| 



1 + 2(|A'|2 + |C/ 



l\2\ 



(|A 



l\2 



\u- 



/|2N2 



2An|[/'| 



and w - Ao - 2\U'\^Ao in the limit \U'\,\A'\ « 1. 
Notice that the bound state has energy below the "local" 
gap magnitude Aq -t- A ~ Ao(l — |t^'P), indicating that 
the formation of in-gap bound state is robust to local gap 
distortion. 
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Impurity averaged single particle density of states 

To have a more quantitative feeling of the effect of im- 
purities we compute the single-particle density of states 
in our model with Uu = U22 = aUi2 = all and A = 0. 
In this case the matrix Green's function G{lu; U) = 
(1 + Go{uj)Mi{U))~^Go{uj) can be evaluated exactly and 
the impurity averaged Green's function 

{G{lo)) = j dUG{LO]U)P{U) (14) 

can be evaluated for given impurity potential distribu- 
tion P{U). Contrary to self-consistent Born type calcu- 
lations the calculation here is valid only in the limit of low 
concentration of impurities where interference effects be- 
tween different impurity scattering events are negligible. 
We find that {Gi2{uj)) = (6*21(0;)) = for even distribu- 
tion P{U) ~ P{—U) and only intra-band Green's func- 
tion survives impurity average. (Inter-band Green's func- 
tion will contribute when evaluating two-particle corre- 
lation functions.) The trace of the total electron Green's 
function is 

TrG(o)) = -2rcj[l + C/'2(l + a^)] 

a{uj)U'^ + b{uj)U'^+c{oj)' ^ ' 

with 

a(w) = {\~o?fSr(uj)S2{u:). 

h{uj) = -2 {u? - IA01A02I cos J - a^Sx{uj)S2{i^)\ , 

c{lo) = Sx{lo)S2{lo). (16) 

The single-particle density of states given by p(lo) = 
3TrG(u;) is evaluated numerically for s± superconduc- 
tors (8 = tt) with P{U) = y|:e-"^" for two differ- 
ent values of a = 0.5 to 3 and a = 0.4 with |Aq2| = 
0.5, 1 a;,! I = 1 and A^i(O) = iV2(0). The results are shown 
in fig. 111). The density of states for a ~ 0.95,1.05 and 
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Figure 1; The averaged density of states p{uj) for a = 0.5, 1, 3, 5 
and a = 0.4 with |AJ,2| = 0.5, |A^^| = 1 and Ni{0) = 7V2(0). 

< u! < A02 is also shown in fig.© for comparison. We 



find that nonzero density of state piw)) is induced inside 
the gap with nonzero spectral weight on the Fermi sur- 
face in general. The precise form of p(w) depends also 
on Uii,U22 and the distribution of impurity scattering 
strength P{U). The in-gap spectral weight increases with 
decreasing a and shifts to lower energy with decreasing a, 
indicating that in-gap states are strengthened by strong 
inter-band scattering but suppressed by intra-band scat- 
tering. Notice also that for a>l,r<l/2 and no 
Lu = bound state exists! The in-gap states have en- 
ergy oj > u!c ^ 0.2 for a = 1.05 as shown in fig. ([2]). We 
cautioned that we haven't included the induced changes 
in gap functions A's in our calculation. Thus we expect 
that the near gap edge behavior in our calculation may 
not be reliable but the deep in-gap behavior should be 
qualitatively correct. 




0) 



Figure 2: The averaged density of states p(tu) for a = I and 
a = 0.95, 1.05 with |A^2| = 0.5, \A'f^^\ = 1 and 7Vi(0) = N2(0). 

Our result indicates that non-magnetic impurities is 
a relevant perturbation to the physics of multi-band su- 
perconductors with frustrated sign between order param- 
eters. They introduce a "ferromagnetic" Josephson cou- 
pling between the two superconducting order parameters 
which disfavors the s± state and introduce in-gap bound 
states in the quasi-particle spectrum. The energy spec- 
trum of the in-gap states depend on the nature of the 
impurity potential and is "non-universal" . 



Origin of the in-gap states 

It is important to understand why in-gap bound states 
exist so easily when the relative phase between the two 
superconducting order parameters is tt, independent of 
the microscopic details of the system. The robustness 
of the in-gap bound states can be understood if we 
notice that the approximate mean-field Green's func- 
tions (j4|) we used in our calculation is the exact Green's 
function of a corresponding one-dimensional supercon- 
ductor problem, if we linearize the fermion spectrum 
e.j: ~ 11 ±Vfi{k — kpi) in action ([1]). In this case the 



5 



density of state Ni{e) ^ dk/de becomes constant and 
Eq. (|4]) becomes exact. In this representation, the impu- 
rity introduces finite tunneUing probabihty between two 
one-dimensional superconductors, one located on the left 
of the impurity, and the other one on the right when 
Ui2 4^ 0. This problem has been studied in Refs.Q 
and |15| . where bound states are found to exist when 
the phase difference between the two superconductors is 

(5 = TT. 

The existence of in-gap bound state can be understood 
by noting that the one-dimensional Bogoliubov equation 
with a linearized electron spectrum around the Fermi sur- 
face is essentially a Dirac equation for spinless fermions 
at one-dimension IJ] . In particular, the gap function A 
becomes the mass term in the Dirac fermion representa- 
tion. Thus the problem is mathematically equivalent to a 
tunnelling problem between two species of Dirac fermions 
with masses of opposite sign. For perfect tunnelling, it is 
known that a = mid-gap state exists if the masses of 
the Dirac fermions have opposite sign at the two sides of 
the tunnelling barrier because of the topological structure 
of the problem [31. The bound states split into two with 
energies ±w > when a tunnelling barrier exists 1^, 17 1 
and eventually merge into the continuum spectrum when 
the tunnelling barrier is high enough. This is in qualita- 
tive agreement with what we find here when = 0. 
In addition, we also observed that the induced A^ is 
small and does not affect the in-gap bound states in the 



limit|C/{2|,|A^ 



<< 1 in the case I A 



on 



A02I and 



Vii = U22 = 0. 

Lastly, we comment that the effect of magnetic impu- 
rity can also be included straightforwardly by introducing 
a magnetic scattering term V^j"^ {r) ^ 5{r)<7VijSaa' in our 
calculation. It is known that magnetic scattering intro- 
duces in-gap bound states even for single-band supercon- 
ductors and this feature remains for two-band supercon- 
ductors without phase frustration. Therefore magnetic 
impurities do not induce qualitative differences between 
two-band s-wave and s±-wave superconductorsfisj which 
is why we concentrate only on non-magnetic impurities 
in the present paper. 



Summarizing, by studying the one-impurity problem 
carefully, we show in this paper that (non-magnetic) 
impurity is a relevant perturbation in two-band s±-wave 
superconductors where the order parameters are of 
opposite signs. The suppression of s± state and gener- 
ation of in-gap bound states in quasi-particle spectrum 
are natural phenomena associated with non-magnetic 
impurities and should be taken into account carefully 
in understanding the electronic properties of Iron-base 
superconductors. 
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